arXiv: 1505.01212v3 [math.AP] 8 Aug 2015 


Stationary solutions of the Vlasov-Fokker-Planck 
equation: existence, characterization and 
phase-transition 


M.H. Duong‘S’*, J. Tugaut'’ 

Mathematics Institute, University of Warwick, Coventry CV4 7AL, UK. 

^ Universite Jean Monnet, Institut Camille Jordan, 23, rue du docteur Paul Michelon, 
CS 82301, 42023 Saint-Etienne Cedex 2, France. 


Abstract 

In this paper, we study the set of stationary solutions of the Vlasov-Fokker- 
Planck (VFP) equation. This equation describes the time evolution of the 
probability distribution of a particle moving under the influence of a double¬ 
well potential, an interaction potential, a friction force and a stochastic force. 
We prove, under suitable assumptions, that the VFP equation does not have 
a unique stationary solution and that there exists a phase transition. Our 
study relies on the recent results by Tugaut and coauthors regarding the 
McKean-Vlasov equation. 

Keywords: 

Invariant measure, Vlasov-Fokker-Planck equation, McKean-Vlasov 
equation, stochastic processes. 

2000 MSG: 60G10, 35Q83, 35Q84. 


1. Introduction 

1.1. The Vlasov-Fokker-Planck equation 

We consider the following Vlasov-Fokker-Planck (VFP) equation, 

dtp = - divg + divp (^piVqV + V q'l) * p + 7^^)) + fkTKpp. (1) 
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In this equation, the spatial domain is with coordinates {q,p) G x 
The unknown is a time-dependent probability measure p: [0,T] —)• P(R^'^). 
Subscripts as in divg and Ap indicate that the differential operators act only 
on those variables. The functions V = V{q) and "0 = 'ip{q) are given. The 
convolution ijj^p is defined by {'ip*p){q) = /jg 2 d '4’{q~q')p{q'^p') dq'dp'. Finally 
7 , k and T are positive constants. 

Equation ([T]) is the forward Kolmogorov equation of the following stochas¬ 
tic differential equation (SDE), 


dQ{t) = P{t) dt, 


dP{t) = -VV{Q{t)) dt-Vip* pt{Q{t)) dt--f — dt+ J2'ykTdW{t), ( 2 ) 

m 

where pt is the law of (Q^, Pt). This SDE models the movement of a particle 
with mass m under a fixed potential V, an interaction potential ip, a friction 
force (the drift term —7 — dt) and a stochastic forcing described by the d- 
dimensional Wiener measures W. In this model, 7 is the friction coefficient, 
k is the Boltzmann constant and T is the absolute temperature. 

Eq. ([T]) and system ([ 2 ]) play an important role in applied sciences in 
particular in statistical mechanics. For instance, they are used as a simplihed 
model for chemical reactions, or as a model for particles inte racting through 
Coulo mb, gravitational, or volume exclusion forces, see e.g., |Kra40 . NPT84 . 


BD95j. Eq. ([T]) (and related models) has been studied intensively in the 
literature by many authors from vario us points of view, see e.g. |Deg 8 fi . 


BD95 . BGMIO . DPZI 3 I. DPZ14 . Duoisj ] and references there in. In particular, 
invariant probab ilities of Eq. ([1]) have been investigated in |Dre87 . BGMld 
(see also Duol5j |b However, in these papers, the potential V is assumed to 
be either bounded or globally Lipschitz or convex. As a result, there is a 
unique stationary solution. In this paper, we show that when the potential 
V is unbounded, non-convex and not globally Lipschitz, of which a double¬ 
well potential is a typical example, non-uniqueness and phase transition can 
occur. Herein, we characterise the set of stationary solutions in such a case. 
Our study relies on the recent results by Tugaut and co-authors about the 
McKean-Vlasov diffusion by showing that the set of stationary solutions of 
the Vlasov-Fokker-Planck equation is related to that of the McKean-Vlasov 
equation. 
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1.2. Normalization 


We first write ([T]) in dimensionless form. The non-dimensionaliz ation for 


Eq. ([T]) has been done previously in the literature, see for instance LRSIO 


Section 2.2.4], For the sake of convenience, we perform it here. By setting 


^ _ mL _ ~ 

q =: Lq, p =: - p, t =: rt 

T 

and 

Viq) ^(^) =• PiPN,t) =: ^^^^j^pip,q,t), 

where L is the characteristic length scale, and r := ^ is the relaxation time 
of the particle dynamics. Then the dimensionless form of the Vlasov-Fokker- 
Planck equation is (after leaving out all the tilde) 


dtp = - divg (^pp^ + divp (p{VqV + Vq^j * p + p)^ + XApp. (3) 

where A = kTT^m~^L~^ is the dimensionless diffusion coefficient. 

In this paper, we are interested in stationary solutions of Eq. ([3]), i.e., 
solutions of the following equation 


K[p](p) = 0, 


( 4 ) 


where 

K[/i](p) := - divg divp [p{VqV + Vqip * /i + p) j + AApp (5) 

for given p G Note that for a given p, the operator K[p]{p) is linear 

in p. This can be seen as a linearised operator of K[p](p). 


1.3. Organisation of the paper 

The rest of the paper is organised as follows. In Section [21 we state our 
assumptions and provide a characterization via an implicit equation for a 
solution of Eq. (jl]). In Section [3] we present main results of the paper which 
prove the existence, (non-) uniqueness and phase transition properties of such 
stationary solutions. 
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2. Characterization of invariant probabilities 

In this section, we characterize solutions of Eq. (|1]). 

First of all, we consider the following assumptions: 

Assumption 1 (Assumptions for the potential V). 

The potential V satishes the following assumptions. 

(VI) V is a smooth function and there exists m eW and C 2 m > 0 such that 
lini|a,|^+oo = C 2 m, where | • | denotes the Euclidean norm. 

(V2) The equation W{x) = 0 admits a hnite number of solutions. We do 
not specify anything about the nature of these critical points. However, 
the local minima where the Hessian is positive will be denoted by Oq. 
(V3) V(x) > Ca\x\‘^ — C 2 |xp for all x G with C 2 , (^4 > 0. 

(V4) lim Hess V(x) = +00 and Hess V(x) >0 for all x ^ K where A is a 

|xK±oo 

compact of which contains all the critical points of V. 


Assumption 2 (Assumptions for the interaction potential 'ip). 

The interaction potential ip satishes the following assumptions. 

{ipl) There exists an even polynomial function G on M such that ip{x) = 
G(|x|). And, deg(G) =: 2n > 2 . 

{ip2) G and G" are convex. 

Ii/j3) G(0) = 0. 

The simplest example (most famous in the literature) is that V (x) = ^^ 
(i.e., V is a double-well potential) and ip{x) = |x^ for some a (i.e., ip is a 
quadratic interaction). 


Proposition 1. Suppose that Assumption 1 and Assumption 2 hold. If 
there exists a solution p^o E C) L°° of Eq. (jl]) then 


1 (p 


Poc{q,p) = Z^ exp 

where Zx is the normalizing constant 

1 


+ V{q) + Ip * Poo{q) 


Zx = 


exp 


A V2 


+ y{q)+'^* Poo{q) 


dqdp. 


( 6 ) 

( 7 ) 


Conversely any measure whose density satishes dH]) is invariant for ([3]). 
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Let us remark that we use the convexity at inhnity of ^|J in order Zx to be 
hnite and poo to be in fl 


Proof. The idea of the proof has appeared in [Dre87l |. where the authors 
study the Vlasov-Fokker-Planck equation but with different scaling and as¬ 
sumptions. The proof is divided into two steps. 

Step 1. We hrst consider the linearised equation 

K(p) ;= K[/x](p) = 0, (8) 

where p G is given. We prove the following assertion: dehne 

u{q,p) := Cf^exp + V (q) + ^ , (9) 

where Cx is the normalisation constant so that UnllLi = 1, and 


A\={v ■. 


D2d 




n(-,p) G Vg G and 


f := V ■ u satisfies / G H^{(1 + \p\ -|- \h\) dqdp), Apf G L^\, 


where h{q) := VqV{q) + Vq'ip * P- Then u is the unique solution in A of the 
linearised equation (|8]). 

Note that under the assumption that V and are sm ooth, t he linearised 
operator is hypo-elliptic, see for instance Hor67l | and DMlfll . Section 1]. 
Hence we know a priori that all solutions of Eq.(IH]) are smooth. Therefore, 
all the derivatives in this proof can be understood in the classical sense. 

We now prove this assertion. By the assumptions on V and fj, V{q) + 
* /i(g) behaves like a polynomial of order fc > 2 of |g| at inhnity. This and 
([9]) imply that 

= / m(1 + \p\ + \h\) dqdp < CX), 


ll^'“niL2((l+|p| + |h|)dqdp) — 


1 


m(|pP -F |hp)(l -h IpI -h \h\) dqdp < oo, 


~ dqdp<oo. 
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Therefore u E A. Since — divg{up) + diVp(M(Vqi^ + Vg'ip * p)) = divp{up) + 
XApU = 0, it follows that K[p]{u) = 0. Now assume that Eq. (|S]) has another 
solution V E A and ||n||Li = 1. Let f := v ■ . We have 

- div,(np) + diVp{v{'\/qV + Vq'0 * fi)) 

= divg(/p) + diVp(/(Vql/ + Vqtl^ * fi))], 

and 

divp (^vp + WpV^ = divp (^vp + XVp{u /) j 

= divp (^vp + X{uVp{u~^^‘^f) + VpU ■ 

= Adivp (uVp{u~^^‘^f)y 

Dehne Q/ := —Then from the above calcula¬ 
tion, we get 

Qf = -[- divgifp) +diVp{f{VqV AVg-ilj * p))] - Xu~^^‘^ divp (uVp{u~^^‘^ f)y 

Therefore, by multiplying by / and integrating over we obtain 
{Qf, /)l 2 = ^ / [divg(p/^) - diVpifiVgV + Vgip * p))] dqdp 
— xf divp (uVp{u~^^‘^f)]f dqdp 

= l [ [diVg(p/^) - diVp{f{VqV + Vgip * p))] dqdp 

^ JK2d 

-xl A\^p{u~^l’^f\vSIp{u~^l’^fyydqdp 
+ X f u(Vp{u~^^‘^f)) dqdp 

jR2d 2 

= x [ u(Vp{u~^^‘^f)') dqdp. 

Note that in the above computations, as we show in Remark [1] that due to 
divergence theorem and the fact that v E A, the hrst three integrals vanish. 

Since Qf = 0, it follows that Vp(M“^/^/) = 0, i.e., u~^l'^f = g{q) for some 
function g. Hence v = vfl‘^ f = u ■ g{q), and 0 = K(n) = —up ■ Vgg{q). It 
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implies that Vqg{q) = 0, i.e., is a constant. Since ||n||Li = 1, we obtain 
that S' = 1, i.e. v = u. In other words, Eq. ([HD has n as a nniqne solntion in 
A and ||m||i = 1. 

Step 2. Snppose that poo G fl L°° is a solntion of Eq. (jlD- Therefore, poo 
solves the eqnation K[poo](^^) = 0. According to Step 1 , this eqnation has a 
nniqne solution given by 


u = exp 




Hence z/ = poo, he., Poo satishes (ED- The reverse assertion is obvious due to 
the convexity at inhnity of V and of '0- hH 


Remark 1. We verify here that the hrst three terms in (Q/, /)l2 vanish. 
We present here for the hrst term only, since the computations for the other 
ones are similar. By dehnition of A, we have 

IIp/^IIli = / \pf\ dqdp < \\f\\l2pi+\p\+\h\)dqdp) < oo, 

jR^d 

and 

II divg(p/^)||ii = 2 /" |/||p-Vq/|dgdp< 2 /" |/| |p| |Vg/| dgdp 

jR^d J]R2d 

< 2 f / \ff\p\ dqdp) ' ( [ \Vqff\p\ dqdp) ' 

KJR'^d J \JR?d J 

^ ‘^\\f\\L^{{l+\p\ + \h\)dqdp)\\'^f\\L'^{{l+\p\ + \h\)dqdp) < OO. 

Therefore, by the divergence theorem, see for instance [GGSinl . Section 4.5.2], 
we obtain divq(p/^) dqdp = 0 , which is the desired equality. 


3. Main results 

In this section, we assume that Assumption 1 and Assumption 2 are 
fulhlled. 

Theorem 1. We consider a measure poo on x M'^. If it is an invariant 
probability for (ED then q J^d Poo{q,p)dp is an invariant probability of 

dX{t) = -VV{X{t)) dt-Vfj* Pt{X{t)) dt + VndW{t), (10) 

where pt is the law of X(t). 
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Proof. Denote by p^o the first marginal of poo, he., Poo{q) = J^d Poo{q,p)dp. 
Suppose that poo is an invariant measure for ([3]). According to Proposition [H 
Poo satishes ([6]), i.e., 


Poc{q,p) 


exp 

1 

A 

+ ^(q) + ^*Poc(q)^ 


Im 2 d exp 


[t + + '0 *Poo(g)) 

dqdp 


1 P_ 
A 2 


L 


.1 p_ 


X 


dp 


exp 

-h 

^V(g) -f'0 *Poo(g)) 


!R<i exp 

1 

A 

[^v(g) + '0*Poo(g)) 

dq 


It follows that 


( 11 ) 


Poo{q)= / pociq,p)dp 


exp 

-{{y{q) + '4^ * Poo{q)) 


exp 

-j{yiq) + '<P*poc{q)) 

dq 


According to HTlOal . Lemma 2.2] poo is a stationary measure of the 
McKean-Vlasov SDE 


dX{t) = -VV{Q{t)) dt-V'f* dt + V^dW{t), (12) 

where pt is the law of X{t). Indeed, the convexity at inhnity of V and of 
provides the sufficient conditions of integrability for poo to be an invariant 
probability. This concludes the proof of this theorem. Note that the forward 
Kolmogorov equation associated to the McKean-Vlasov SDE is given by 


dtpt = div[pt(VV + pt)] + AApt. (13) 

□ 


Theorem [T] establishes a one-to-one correspondence between invariant 
measures of the McKean-Vlasov equation flTdl) and that of the Vlasov-Fokker- 
Planck equation ([3]). We therefore can deduce various results on the exis¬ 
tence, (non)-uniqueness and characterization for the latter from the known 
results obtained by Tugaut and co-authors for the former. 

The hrst result corresponds to the existence of an invariant probability 
under Assumptions 1 and 2. 





















Proposition 2. For any A > 0, there exists a stationary measure - that is a 
solution of Eq. (|1]) . 


This is a consequence of Proposition 3.1 in |Tugl4bl |. The idea is to consider 


the free-energy functional. By using some compactness argument, we are able 
to prove that a subsequence of the law of X{t) admits an adherence value 
which is an invariant probability, where X is the McKean-Vlasov diffusion. 
Then, it proves the existence of such an invariant probability. 

Next results do describe the invariant probabilities. 

First one deals with the case in which both V and "0 are even functions. 
In this case, by using Schauder hxed point theorem, we are able to obtain 
the existence of an invariant probability whose density with respect to the 
Lebesgue measure is even. 


Theorem 2. If both V and are even, there exists a symmetric invari¬ 
ant probability (that is a measure whose density with respect to the Lebesgue 
measure is even). 


This is a consequence of Theorem 4.5 in HTlOa 


Next proposition establishes that there is an invariant probability around 
each well, providing an additive condition on the potentials. 


Proposition 3. Here, d = 1. We assume that the interacting potential 
ip is quadratic: := Let oq be a critical point of V such that 

a -f- V'\af) > 0 and 


a > 2 sup 

x^ao 


V(ao) - V(x) 
(ao - xf 


(14) 


Thus, for all d g] 0 ; 1[, there exists Aq > 0 such that for all A < Aq, Diffusion 
(|3]) admits an invariant probability p^o satisfying 



qpoo{q,p)dqdp- ao + 


V^^\ao) 


4V"(ao) (a + K"(«o)) 


A 


< 6X 


This is a consequence of Proposition 1.2 in Tugl4al . 

Let us explain briefly what Condition flTT)) means. When the temperature 
is small, each stationary measure which is not symmetric concentrates to 
a Dirac measure around a critical point. However, the well needs to be 
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sufficiently deep in the following sense: for any x 7 ^ Oq, one needs that 
l/(x) + I (x — Oo)^ > V^(ao)) which is exactly equivalent to Condition (IT^ . 

In the following theorem, we take assumptions close to the ones in Propo- 
sition|3l Indeed, the interacting potential -0 is assumed to be quadratic. How¬ 
ever, by using some convexity assumptions on the derivatives of V, we will 
show that there are exactly one or three invariant probabilities. Moreover, 
we can simulate if we are in the uniqueness or in the thirdness case. 

Theorem 3. Here, d = 1. We assume that 

V{x) = -f V with deg(H) =: 2q . (15) 

2 ^ (2P)! 

And, ip{x) := |x^. Thus, there exists Ac > 0 such that: 

• For all X > Ac, Diffusion ([3]) admits a unique invariant probability, 
which is symmetric. 

• For all A < Ac, Diffusion ([3]) admits exactly three invariant probabili¬ 
ties. 


Moreover, Ac is the unique solution of the equation, where z is the unknown, 

2zP~^\v^^P'>(0)\ 


2a 




(2p)\ 


l2'2Py2p 


dy = 0. 


(16) 


This is a consequence of Theorem 2.1 in |Tugl4a |. The global idea is the 
following. The study of the invariant probabilities is equivalent to the study 
of the zeros of a function from M to M. The behaviour of the function does 
depend on the value A. If A > Ac, this function has a unique zero and if 
A < Ac, the function is increasing then decreasing on M+ so it admits three 
zeros. 

A similar idea is used to obtain the following result. 


Proposition 4. Here, d = 1. We assume that if is quadratic: ijj{x) := |x^. 
Thus, for any a > 0, there exists a critical value Ao(a) such that Diffusion 
([3]) admits a unique invariant probability provided that A > Ao(a). 


This is a consequence of Proposition 2.4 in |Tugl4a |. In Proposition 


m we provide a setting such that there is the uniqueness of the invariant 
probability. The main difference with Theorem [3] is that we do not assume 
that the derivative of V is convex. 
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Theorem 4. Let Qq be a point where V admits a local minimum such that 


V{x) +'ifjix — Oq) > V {oq) for all x ^ Oq . (17) 

Then, for all k > 0 small enough, there exists Aq > 0 such that VA g]0; Ao[, 
the diffusion (E]) admits a stationary measure poo satisfying 



\\q - poc,{q,p)dqdp < K 


2n 


This is a consequence of Theorem 2.3 in [Tugl4b|. The global idea is the 


same than the one of Proposition E] but the technical material is different. 
Indeed, we use the free-energy funct ional. 


Mo re generally, all the result s in jHTlOal. iHTlObl. lHT12l. IMcKGGI. IMcKGTI . 


TuglOl. iTuglll. lTugl4ai lTugl4bl | hold. 
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